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Abstract 

It is proposed that convection driven dynamos operating in planetary cores could be 
oscillatory even when the oscillations are not directly noticeable from the outside. 
Examples of dynamo simulations are pointed out that exhibit oscillations in the 
structure of the azimuthally averaged toroidal magnetic flux while the mean poloidal 
field shows only variations in its amplitude. In the case of the geomagnetic field, 
global excursions may be associated with these oscillations. Long period dynamo 
simulations indicate that the oscillations may cause reversals once in a while. No 
special attempt has been made to use most realistic parameter values. Nevertheless 
some similarities between the simulations and the paleomagnetic record can be 
pointed out. 
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1 Introduction 



The origin of geomagnetic reversals is a much debated subject among scientists 
in the fields of paleomagnetism and dynamo theory. There is general agreement 
that a detailed understanding of reversals is a key issue of geodynamo theory. 
In this connection also the problem of global excursions of the geomagnetic 
field in which the dipole strength reaches te mporarily unusually l ow values 
has been discussed and it has been suggested Woell & Cox I (1 1 9 72l ). I ffoffm an 



Email addresses: busse@uni-bayreuth.de (F. H. Busse), 
r.simitev@niaths.gla.ac.uk (R. D. Simitev). 



Phys. Earth Planet. Inter., 168(3-4), 237, 2008, doi:10.1016/j.pepi.2008.06.007 



(119811) ■ see als o IMerrill et al. I (119961 ) and later papers by iGubbins I (119991 ) and 
Wicht I tOO^ ) that global excursions are aborted reversals. Not all recorded 
excursions are global one s, but global excursion s still appear to occur more 
frequently than reversals, wanaereis et al. I (119971 ) identified at least six global 
excursions in the last about 800 ky since the Brunhes/Matuyama reversal, but 
in the last d e cade many more global excursions have been found according to 
(120061 ). In this letter we wish to support the notion that global 



Lund et al. 



excursions and reversals originate from the same mechanism. An oscillatory 
dynamo process that manifests itself primarily in the toroidal component of 
the magnetic field will be proposed as such a mechanism. Indeed, from the 
perspective of the oscillations, excursions must be considered as the normal 
behavior, while a reversal represents an exceptional excursion in which the 
mean poloidal field is perturbed more strongly that it can recover from its 
low-amplitude state only with the opposite sign. 

Traditionally the geodynamo is regarded as a stationary dynamo in contrast to 
the solar dynamo which exhibits a 22-year period. Dynamo simulations have 
shown, however, that in rapidly rotating spherical fiuid shells with significant 
differential rotation often oscillatory dynamos are found. That dynamo oscil- 
lations may not be visible from the exterior of the conducting fluid sphere 
has been pointed out previously ^Busse & Simitev I (120061 )). The present let- 



ter intends to demonstrate how oscillations can lead to global excursions and 
more rarely to reversals. While the simulations are based on the fundamental 
equations governing the generation of magnetic fields by convection fiows in 
rotating spherical shells, only a minimum of physical parameters is introduced 
and a faithful modeling of the Earth's core has not been the primary goal. 



2 Mathematical formulation 



We consider a spherical fiuid shell of thickness d rotating with a constant 

angular velocity Q. It is assumed that a static state exists with the temperature 

distribution Ts = Tq — (3d'^r'^/2. Here rd is the length of the position vector, 

r, with respect to the center of the sphere. The gravity field is g = ~d-fr. 

In addition to the length d, the time d'^/u, the temperature u'^/'yad'^ and 

the magnetic fiux density u^hqY^'^ /d are used as scales for the dimensionless 

description of the problem where u denotes the kinematic viscosity of the fiuid, 

K its thermal diffusivity, g its density, a its coefficient of thermal expansion 

and /i is its magnetic permeability. The Boussinesq approximation is assumed. 

Accordingly, the velocity field u as well as the magnetic fiux density B are 

solenoidal vector fields for which the general representation in terms of poloidal 

and toroidal components can be used, 

u = V X (Vw X r) + Vw X r , (la) 

B = V X (V/i X r) + S7g xr . (lb) 
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By multiplying the (curl)^ and the curl of the equation of motion and of the 
induction equation by r, we obtain four equations for v and w and for h and 
g. These four equations together with the heat equation for the dimensionless 
deviation 6 from the static temperature distribution and with the appropriate 
boundary conditions represent the basis for the mathematical description of 
the evolution in time of thermal convection in the rotating spherical shell and 
of the magnetic fi eld generated by it . Sinc e these equations ha v e been given in 
previous papers \Simitev & Busse I feoosh . \Busse & Simitev \ fl2006h ). we list 



here only the dimensionless parameters, the Rayleigh number R, the Coriolis 
number r, the Prandtl number P and the magnetic Prandtl number P^, 

i^^^, r^?^. P^^, (2, 

UK U K A 

where A is the magnetic diffusivity. We assume stress-free boundaries with 
fixed temperatures and use the radius ratio ri/vo = 0.4, 
V = dl^v = dr{w/r) = 6 = 

at r = Tj = 2/3 and r = To = 5/3. (3) 

For the magnetic field an electrically insulating outer boundary is assumed 
such that the poloidal function h must be matched to the function h^^^ which 
describes the potential field outside the fluid shell 

g = h- h^'^ = dr{h - h^^^) = QXr = ro = 5/3. (4) 

In order to avoid the computation of h and g in the inner core, r < r^, we 
assume either an electrically insulating inner boundary, 

g = h- h'^'^ = dr{h - h^^^) = at r = ri = 2/3, (5) 

or a perfectly conducting inner core in which case the conditions 

h = dr{rg) = at r = ri = 2/3 (6) 

must be applied. The numerical integration of the equations for v, w, 0, h and 
g together with boundary conditions (3), (4) and ( 5) or (6) proceeds with 



the pseudo-spectral method as described hy \Tilgner \ (119991 ) which is based on 
an expansion of all dependent variables in spherical harmonics for the 6, ip- 
dependences, i.e. h = Y, H^r, t)Pr{cos 9) exp{imip} (7) 

l,m 

and analogous expressions for the other variables, v, w, and g. P{^ denotes 
the associated Legendre functions. For the r-dependence expansions in Cheby- 
chev polynomials are used. Azimuthally averaged components of the fields 
V, w, 0, h and g will be indicated by an overbar. For most computations to 
be reported here a minimum of 33 collocation points in the radial direction 
and spherical harmonics up to the order 96 have been used. But this high 
resolution was not needed in all cases. Instead of the time t based on the vis- 
cous time scale we shall use in the following the time t* = t/ Pm based on the 
magnetic diffusion time, cP/X. 
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3 Oscillations of the Toroidal Magnetic Flux 



Even in their turbulent state of motion, convection flows outside the tangent 
cylinder which touches the inner core boundary at its equator remain essen- 
tially symmetric with respect to equatorial plane as is evident from flgure 1. 
For this reason dynamo solutions characterized by an axial dipole correspond 
to a mean azimuthal magnetic flux that is antisymmetric with respect to the 
equatorial plane. Oscillations of these axisymmetric flux tubes originate from 
the creation of a pair of new flux tubes with opposite signs at the equatorial 
plane which grow and push the older flux towards higher latitudes as shown in 
flgure 2. This process is strongly dependent on the differential rotation which 
is prograde at larger distances from the axis and retrograde at smal ler ones. 



The o scillations can be desc ribed by Parker's dynam o wave model 1] Parker 



( 1l955l )) as has been done by iBusse & Simitev I ( l2006l ). In the present case of 



flgure 2 the oscillation is modifled in two respects. First, the mean toroidal 
fleld becomes nearly quadrupolar, i.e. symmetric about the equatorial plane, as 
the amplitudes of the mean poloidal fleld and of the differential rotation reach 
their minimum values. Secondly and more importantly, the mean poloidal fleld 
participates in the oscillation only as far as its amplitude varies. In the case 
of flgure 2 its amplitude decays and reaches a minimum around t* ^ 1.6 at 
which time a magnetic eddy emerges with the opposite sign of the given mean 
poloidal fleld. Usually this eddy drifts outward and dissipates as it reaches the 
surface of the conducting region such that the original poloidal fleld prevails. 
Now a relatively long time passes before the process repeats itself and new 
toroidal flux emerges at the equatorial plane. In contrast to the thinner flux 
tubes of dynamos at higher Prandtl numbers which exhibit a more sinusoidal 
oscillation as shown in section (b) of flgure 3, the oscillation in the present 
case resembles more a relaxation oscillation as shown in section (a) of flgure 
3. The amplitudes Hi, in this flgure are assumed at the mid-radius of the 
fluid shell, but usually does not differ much from the dipole component 
describing the magnetic fleld outside the fluid shell. 

While the process visualized in flgure 2 shares several features with global 
excursions, it may also give rise to reversals. These happen in some cases 
when the emerging eddy with the opposite sign of the poloidal fleld replaces 
the latter as shown in flgure 4. This situation occurs most likely if the eddy 
with the opposite sign emerges near the equatorial plane such that it splits 
the original fleld into two parts. It is remarkable that the reversed poloidal 
fleld appears flrst at l ow latitud e s as h as also been observed in the case of 



geomagnetic reversals n Clement I (120041 )). Note that the radius r = + 1.3 



corresponds approximately to the Earth's surface. The occurrence of a reversal 
seems to be promoted by a particularly strong equatorially symmetric toroidal 
flux as appears to be indicated by the correlation between reversals and relative 
high absolute values of the coefficient in sections (b) and (c) of flgure 3. We 
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note in passing that Li et al. (2002), propose a reversal mechanism in which 
the quadrupole mode grows, exceeding the dipole mode before the reversal in a 
manner similar to what happens near the minimum of the oscillations shown in 
figure 2. In contrast, however, our dynamo solutions do not alternate between 
high- and low-energy states, nor do they exhibit a broken columnar vortex 
structure of the velocity field. 

The examples discussed so far all correspond to a single set of parameter val- 
ues. In particular condition (6) for a highly electrically conducting core has 
been used. In order to demonstrate the robust nature of the mechanism of 
global excursions and reversals, we show in figure 5 a sequence of plots ex- 
hibiting a reversal from a dynamo simulation with a quite different set of 
parameters for which condition (5) instead of (6) has been applied. The oscil- 
lations occur somewhat less regularly in this case as is evident from the time 
series of the amplitude of the axial dipole component shown in section (c) of 
figure 3, but the average period is again close to half a magnetic diffusion time. 
A common property of the oscillations is that the quadrupolar components of 
the axisymmetric magnetic field play a significant role. In this respect some 
similarity may be noted with the oscillations displayed in figures 12 and 13 of 



Busse & Simitev I (120061 ). 



Although the inner core does not participate in the oscillations in either of the 
boundary conditions (5) and (6), we expect that the use of a vanishing jump 
of the electrical conductivity at r = will not affect t he resul t s sign ificantly. 



As has been observed i n the dynamo simulations of I Wicht I (120021 ) and of 



Simitev & Busse I (120051 ). because of its small volume the inner core does not 



appear to have a significant effect on the dynamo process. 



4 Discussion 



In selecting the dynamo cases displayed in figures 2, 4 and 5 we have em- 
phasized a high value of r and a reasonably high value of R for which the 
available computer capacity allows to obtain time records extending over 
many magnetic diffusion times. The critical values of the Rayleigh numbers for 
r = 3 X 10^ and r = 10^ are = 2.35 x 10^ with = 10 and R^ = 1.05 x 10^ 
with rric = 11, respectively. Hence the Rayleigh numbers used for the cases of 
figures 2, 4 and 5 exceed their critical values by nearly a factor of four. The 
corresponding average Nusselt numbers at the inner boundary are Nui = 1.58 
and Nui = 1.73 and the corresponding magnetic Reynolds numbers, defined 
by Rm = PmV^E, are Rm = 210 and Rm = 156, respectively. The Prandtl 
number P = 0.1 was used in both cases since it appears to be a reasonable 
comp romise b e tween the molecular value P = 0.05 estimated for the outer 
core ( Poirier ( 19881 )) and a value of the order unity usually assumed for a 
highly turbulent fluid. Moreover, the choice of a low value of P has allowed 
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us to choose a desirable relatively low value of P„ 



The succe ssful application of Parke r's kinematic model for dynamo waves em- 
ployed by iBusse & Simitev I (120061 ) suggests that the oscillations depend pri- 
marily on the differential rotation and the mean helicity of convection which 
are assumed as given. The modified oscillation considered in the present paper 
is characterized by an extended phase of a dominant equatorially symmetric 
(quadrupolar) mean toroidal field which is responsible for the property that 
the period becomes comparable to the magnetic diffusion time. The variations 
of the amplitude of convection and of the differential rotation seem to be of 
lesser importance. 



Using the depth d ^ 2200 km of the liquid outer core and a typical and 
often quoted value A ~ 2 m^/s we find 0.8 x 10^ years as the magnetic dif- 
fusion time of the Earth's core which corresponds to t* = 1 in the figures 
of this paper. The oscillation period T* ^ 0.5 obtained in the time series 
of figure 3(a) thus roughly equals about 40 ky in the Earth's core. This 
period is quite comparable to the broad maximum in the region of 30-50 
ky that seems to char acterize the spectrurn of th e amplitude varia t ions o f 



the geomagnetic fi eld ^Tauxe & Shackleton I (119941 ). \ Tauxe & Hartl\ ( 119971 ) 



Guyodo & Valet 1 (Il999h) throu g hout the last million years. A more recent 
analysis ( Constable & Johnson I ( ]2005l )) has shed some doubts on the exis- 
tence of such a spectral peak, but still confirms a sharp decrease of the spectral 
power for periods shorter than about 30 ky. We also like to draw attention to 
th e property that the typical separation between global excursions in table 1 
oi \Lund et al. I ( 120061 ) varies between 30 and 50 ky. 



From the reversals exhibited in figures 3, 4 and 5 it appears that the amplitude 
increases more sharply after the reversal than it decays towards the reversal. To 
demonstrate this effect more clearly we have plotted in figure 6 the coefficient 

in proximity of the reversal as a function of time for each of the last 4 
reversals that have been obtained in the cases a) and c) of figure 3. Although 
the asymmetry between the dipole strengths before and after the r eversal is not 
as stro r ig as has been found in th e case of paleomagnetic reversals (j Valet et al. 
^200^ . IGuyodo & Valet I ( l2006l )). a similar effect seems to exist. Since the time 
records of figure 3 do not exhibit this effect very well we have plotted in figure 
6 values of at r = + 0.5 for shorter time periods. In the case R = 850000 

aX r = To is also shown (by solid lines) since it represents the axial dipole 
strength of the potential field outside the fiuid shell. Apart from a small shift 
in time the value of does not vary much as function of the radius. In the 
continuing investigation of the dynamo oscillations it will be attempted to find 
even closer correspondences with paleomagnetic observations. 



The possibility of toroidal fiux oscillations as origin of global excursions and 
reversals proposed in this paper differs from all other mechanisms proposed 
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in the literature for reversals and excursions and resembles more the mecha- 
nisms considered for the solar cycle. In the latter the mean poloidal field fully 
participates, of course, similarly as in the dipole oscillation of figure 10 of 



Busse & Simitev I (120061 ) except for the property that the solar dynamo wave 



propagates towards lower instead of higher latitudes. A comparison of different 
mechanisms for geomagnetic reversals would go beyond the scope of present 
paper and should be postponed until more detailed computational results for 
a wider range of parameters become available. 
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Fig. 1. (color online). Typical structures of the velocity field in the case P = 0.1, 
r = 10^, i? = 4 X 10^, Pm = 0.5 with a perfectly electrically conducting inner 
core. The left plot shows lines of constant u^p in the left half and streamlines 
rsmOdev =const. in the right half, all in the meridional plane. The right plot shows 
lines of constant Ur at r = rj + 0.5 at the time t* = 1.486. Positive and negative 
values are indicated by solid (red online) and dashed (blue online) lines. 
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Fig. 2. (color online). Dynamo oscillation in the case 

p = 0.1, T = lo^ R = 4x lo^ 

Pm = 0.5 with perfectly conducting inner core. The half circles show lines of constant 
Uip. The full circles show meridional isolines of B^p (left half) and of r sin Odgh (right 
half) at times t* = 1.490, 1.538, 1.586, 1.634, (first column, from top to bottom) 
and t* = 1.682, 1.810, 1.954, 2.034 (second column). The times t* refer to figure 
3(a). 
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Fig. 3. (color online). Selected coefficients at r = + 0.5 in the cases (a) P = 0.1, 
r = 10^, i? = 4 X 10^, Pm = 0.5 with perfectly conducting inner core; (b) P = 5, 
T = 5000, R = 600000, Pm = 10 with electrically insulating inner core; (c) P = 0.1, 
T = 3 X 10'^, R = 850000, P^ = 1 with insulating inner core. The coefficient of 
the axial dipole component Hf (axial quadrupole component H2) is indicated by a 
solid/red online (dashed/blue online) line. The coefficient G2 in (a) is indicated by 
a dashed /green online line. 
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Fig. 4. (color online). Magnetic field polarity reversal in the case P = 0.1, r = 10^, 
i? = 4 X 10^, Pm = 0.5 with perfectly conducting inner core. The left column shows 
meridional isolines of (left half) and of r sin 9dgh (right half). The right column 
shows lines Br = const, at r = ro+1.3. The interval between the plots is At* = 0.048 
with the first plot at t* = 0.994 (see figure 3(a)). 
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Fig. 5. (color online). Same as figure 4, but for P = 0.1, r = 3 x lO'^, R = 850000, 
Pm = 1 with insulating inner core. The interval between the plots is At* = 0.07 
with the first plot at t* = 26.155 (see figure 3(c)). 
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Fig. 6. (color online). Time-series of the coefficients of the axial dipole component 
Hf at r = + 0.5 across the last 4 reversals in the cases P = 0.1, r = 10^, 
R = A X 10^, Pm = 0.5 with perfectly conducting inner core (top) and P = 0.1, 
r = 3 X 10^, i? = 850000, Pm = I with insulating inner core (bottom). For the sake 
of comparison, the time series have been translated along the time axis so that the 
polarity transitions occur at t = and —H^ is plotted for every second reversal. In 
both panels, black, red, blue and green color correspond to reversals 1(2) to 4(5) of 
the respective cases in figure 3. In the bottom panel, Hf at r = Tq has been included 
in order to represent the axial dipole strength of the potential field outside the fluid 
sheh. Hf at r = ri + 0.5 ( given by dashed lines) precedes it by about At* ~ 0.04 
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